Electronic properties of a biased graphene bilayer 
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We study, within the tight-binding approximation, the electronic properties of a graphene bilayer 
in the presence of an external electric field applied perpendicular to the system - biased bilayer. 
The effect of the perpendicular electric field is included through a parallel plate capacitor model, 
with screening correction at the Hartree level. The full tight-binding description is compared with 
its 4-band and 2-band continuum approximations, and the 4-band model is shown to be always 
a suitable approximation for the conditions realized in experiments. The model is applied to real 
biased bilayer devices, either made out of SiC or exfoliated graphene, and good agreement with 
experimental results is found, indicating that the model is capturing the key ingredients, and that 
a finite gap is effectively being controlled externally. Analysis of experimental results regarding 
the electrical noise and cyclotron resonance further suggests that the model can be seen as a good 
starting point to understand the electronic properties of graphene bilayer. Also, we study the effect 
of electron-hole asymmetry terms, as the second-nearest-neighbor hopping energies t' (in-plane) and 
74 (inter-layer), and the on-site energy A. 

PACS numbers: 73.20.At, 73.21.Ac, 73.43.-f, 81.05.Uw 
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I. INTRODUCTION 



The double layer graphene system - the so-called bi- 
layer graphene (BLG) - is now a subject of considerable 
interest due to its unusual properties,— — dissimilar in 
large extent to those of the single layer graphene (SLG)»^ 
The integer quantum Hall effect (QHE) is a paradigmatic 
case; characterized by the absence of a plateau at the 
Dirac pointii^ and thus still anomalous, it is associated 
with massive Dirac fermions and two zero energy modesji 

One of the most remarkable properties of BLG is the 
ability to open a gap in the spectrum by electric field 
effect - biased BLGii This has been shown both ex- 
perimentally and theoretically, providing the first semi- 
conductor with externally tunable gapi^^— In the ab- 
sence of external perpendicular electric field - unbiased 
BLG - the system is characterized by four bands, two 
of them touching each other parabolically at zero en- 
ergy, and giving rise to the massive Dirac fermions men- 
tioned above, and other two separated by an energy ±ij^. 
Hence, an unbiased BLG is a two-dimensional zero-gap 
semiconductor i^ii^ At the neutrality point the conduc- 
tivity shows a minimum of the order of the conductance 
quantum 



6.18-25 



a property shared with SLGi^S This pre- 
vents standard device applications where the presence 
of a finite gap producing high on-off current ratios is 
of paramount importance. The fact that a simple per- 
pendicular electric field is enough to open a gap, and 
even more remarkable, to control its size, clearly demon- 
strates the potential of this system for carbon-based 
electronics <22i2^ 



The biased BLG reveals interesting properties on its 
own. The gap has shown to be robust in the presence 
of disorder;^"— induced either by impurities or dilution, 
but is completely absent in rotated (non ^i3-stacked) bi- 
layers, where the SLG linear dispersion is recovered i^2i2^ 
The band structure near the gap shows a "Mexican-hat" 
like behavior, with a low doping Fermi surface which is a 
ringF^i Such a topologically nontrivial Fermi surface leads 
to an enhancement of electron-electron interactions, and 
to a ferromagnetic instability at low enough density of 
carriers i^^i^ In the presence of a perpendicular magnetic 
field, the biased BLG shows cyclotron mass renormal- 
ization and an extra plateau at zero Hall conductivity, 
signaling the presence of a sizable gap at the neutral- 



ity point 



10.12.36 



Gaps can also be induced in stacks with 



more than two layers as long as the stacking order is 
of the rhombohedral-typefiiiii although screening effects 
may become important in doped systems with increasing 
number of layersJ^ Recently, a ferromagnetic proximity 
effect was proposed as a different mechanism which can 
also open a gap in the spectrum of the BLG, leading 
to a sizable magnetoresistive effect)^ Strain applied to 
the biased BLG has also shown to produce further gap 
modulation!^ 

In this paper the electronic properties of a biased BLG 
are studied within a full tight-binding model, which en- 
ables the analysis of the whole bandwidth, validating pre- 
vious results obtained using low-energy effective models. 
The screening of the applied perpendicular electric field is 
obtained within a self-consistent Hartree approach, and 
a comparison with experiments is provided. The effect of 
the bias in the cyclotron mass and cyclotron resonance is 



addressed, and the results are shown to agree well with 
experimental measurements. 

The paper is organized as follows: in Sec. |ll] the lat- 
tice structure of BLG and the tight-binding Hamiltonian 
are presented; bulk electronic properties are discussed in 
Sec IIIH with particular emphasis on the screening cor- 
rection; the effect of a perpendicular magnetic field is 
studied in Sec. IIV[ Sec. |V] contains our conclusions. We 
have also included three appendices: Appendix |X] pro- 
vides details on the calculation of the density asymmetry 
between layers for a finite bias; in Appendix |B] we give 
the analytical expression for the biased BLG density of 
states, valid over the entire energy spectrum; analytical 
expressions for the cyclotron mass obtained within the 
full tight-binding model are given in Appendix [C] 



II. MODEL 

Here we consider only ^i?-Bernal stacking, where the 
top layer has its A sublattice on top of sublattice B of 
the bottom layer. We use indices 1 and 2 to label the 
top and bottom layer, respectively. The unit cell of a 
bilayer has twice the number of atoms of a single layer. 
The basis vectors may be written as ai = a e^; and sl2 = 
a{ex — Viey)/2, where a = 2.46 A. 

In the tight-binding approximation, the in-plane hop- 
ping energy, t, and the inter-layer hopping energy, t^, 
define the most relevant energy scales. The simplest 
tight-binging Hamiltonian describing non-interacting 
TT— electrons in BLG reads: 

2 

Htb = ^i/.+t^^ [aI,,(R)fe2,.(R)+li.c.]+Hy, (1) 

i=l R.cr 

with the SLG Hamiltonian 

H, = -tY^ [a[^(R)6,^,(R) + a[^(R)6,,,(R - ai) 



73/*^^ 


. /^41.48.49 


A/A^-^ 


f'/A^ 




0.03-0.1 


0.04 - 0.07 


0.005 - 0.008 


-0.04 



Table II: Approximate parameter values as obtained in recent 
experiments (except for t' quoted from DFT calculations). 



using the well established Slonczewski-Weiss-McClure 
(SWM) parametrization model^i^ to fit the bands near 
the Fermi energy. The SWM model assumes extra pa- 
rameters that can also be incorporated in a tight-binding 
model for BLG. Namely, the inter-layer second- NN hop- 
pings 73 and 74, where 73 connects different sublattices 
{B1-A2) and 74 equal sublattices {A\~A2 and B1-B2). 
Additionally, there is an on-site energy A reflecting the 
inequivalence between sublattices A1,B2 and B1,A2 - 
the former project exactly on top of each other while the 
latter lay on the hexagon center of the other layer. The 
consequences of these extra terms for the band struc- 
ture obtained from Eq. ([T]) are well known: 73 induces 
trigonal warping and both 74 and A give rise to electron- 
hole asymmetry 1^ 1 "^^1^^ The in-plane second-NN hopping 
energy t' is not considered in the usual tight-binding 
parametrization of the SWM model. Nevertheless, this 
term can have important consequences since it breaks 
particle-hole symmetry but does not modify the Dirac 
spectrum. Typical values are given in Table |lT] as ob- 
tained in recent experiments, except for t' quoted from 
density functional theory (DFT) calculations. 

We are interested in the properties of BLG in the pres- 
ence of a perpendicular electric field - the biased BLG. 
The effect of the induced energy difference between lay- 
ers, parametrized by V , may be accounted for by adding 
Hv to Eq. IT]) , with Hy given by 

Hv = ^ I] [riAiiR] + nsiCR) - 7iA2(R) - nB2{R)] , 



R,r 



R,cr 



a'- (R)6 (R — ao) + h c 1 (2) where n^i(R) and nBiCR) are number operators. 



(3) 



where ai^o-(R) [^i.o-lR)] is the annihilation operator for 
electrons at position R in sublattice Ai (Bi), i — 1,2, 
and spin a. The in-plane hopping t can be inferred 
from the Fermi velocity in graphene vp — tah~^y3/2 w 
10^ms~^)2S yielding t sa 3.1 eV, in good agreement with 
what is found experimentally for graphite4S This value 
also agrees with a recent Raman scattering study of 
the electronic structure of BLG.~ As regards the inter- 
layer hopping tj_, angle-resolved photoemission spec- 
troscopy (ARPES) measurements in epitaxial BLG give 
t± « 0.43 eV,— and Raman scattering for BLG ob- 
tained by micro-mechanical cleavage of graphite yields 
t± « 0.30 eV4i The experimental value for bulk graphite 
is t± « 0.39 eV)^ which means that for practical pur- 
poses we can always assume t±/t ^ 0.1 <C 1. This 
values for t and t± compare fairly well with what is 
obtained from fist-principles calculations for graphite^ 



III. BULK ELECTRONIC PROPERTIES 

Introducing the Fourier components Ci^cr.k and 
^i.cr.k of operators ai^a(R-) and bi^^CR), respectively, 
with the layer index i — 1,2, we can rewrite 
Eq. H]) as iJ = Ek,^ ^'ik-f^k^^,k, where i/jI^^ = 

["i.cr,k' ^i.cr.k' "2.cr,k' ^2.cr,k] ^^ ^ f°™ compoueut spluor, 
and iJk is given by 



-ffk = 



72 -isk 





-ti_ 


tsl V/2 











-v/2 


-isk 


-tj_ 


-tsl 


-V/2 



(4) 



The factor Sk = 1 + e^^'^^ + e^^'^^ determines the matrix 
elements for the SLG Hamiltonian in reciprocal space 



(i_L = 0, F = 0), from which the SLG dispersion is ob- 
tained, eji = ±i|sk|. The resultant conduction (+) and 
valence (— ) bands touch each other in a conical way at 
the corners of the first Brillouin zone (BZ), the K and K' 
pointsi^ This touching occurs at zero energy, the Fermi 
energy for undoped graphene. The 4-band continuum 
approximation for Eq. Q, valid at energy scales E <^t, 
may be obtained by introducing the small wave vector q 
which measures the difference between k and the corners 
of the BZ. Linearizing the factor Sk around the K points 
Eq. (g]) reads 



H 



K 



V 



V/2 


VFpe *'^p 





-ti_ \ 


UFpe*'^p 


V/2 














-v/2 


VYpe'^"^^ 


-ti_ 





VFpe^'^'p 


-V/2 j 



(5) 



where p = ?iq and Lp-p = tan ^{py/px). Around the K' 
points Eq. ([5]) with complex conjugate matrix elements 
defines Hk"^^ 

Equation ([5]) can be further simplified if one assumes 
vfp, V <S^ t±. By eliminating high energy states pertur- 
batively we can write a two-band effective Hamiltonian 
describing low-energy states whose electronic amplitude 
is mostly localized on Bl and A2 sites. Near the K points 
the resulting Hamiltonian may be written as 



^eff - - 



-~V/2 
e-^'^'P^vlp^/ti 



V/2 



(6) 



whereas the complex conjugate matrix elements should 
be taken for a low-energy description around the K' 
points. The two-component wave functions have the form 

In the following we discuss the electronic structure re- 
sulting from the tight-binding Hamiltonian Q , and com- 
ment on the approximations given above by Eqs. ([5|) 
and ©. 
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Figure 1: (Color online) (a) Solution of Eq. ^ ior V = 
t±/2,2t±,4:t±. (b) Ag vs V for various t± values. Energy 
is given in units of t and momentum in units of a~^ . 



It can be readily shown that the minimum of sub-band 
(V) [or equivalently, the maximum of E^~{V)\ oc- 
curs for all k's satisfying 



K 



aiV,t^), 



(8) 



with a{V, ij.) = (y4/4 + t\V^/2)/{V^ + t\) - note that 
dE^^ /de]^ = at the desired extrema. Equation ([5]) has 
solutions for ^/a < 3i (3i is half of the single layer band- 
width). When y/a > 3t the minimum of E^~ (V) occurs 
at the F point. Figure [TJ a) shows the solution of Eq. ([5]) 
around the K point for V = t±/2,2t±,4:t± (around the 
K' point the figure is rotated by tt/B). At low doping 
the Fermi sea acquires a line shape given by the solution 
of Eq. ([H]) , the line width being determined by the dop- 
ing level. As can be seen in Fig. [ija), when V < t± the 
Fermi sea approaches a ring, the Fermi ring, centered at 
the BZ corners. J^^ As V is increased there is an appar- 
ent trigonal distortion showing up, which originates from 
the single layer dispersion in Eq. ([8|). 

The existence of a Fermi ring is easily understood using 
the continuum version of Eq. ([T]), i.e., the eigenvalues of 
Eq. ([5|). This amounts to substituting the single layer 
dispersion in Eq. ([7|) by vpp, which immediately implies 
cylindrical symmetry around K and K' . If we further 
assume that vpp -^ V -^ t± holds, Eq. ([7]) is then well 
approximated by the "Mexican hat" dispersion^ii 



Let us briefiy discuss the electronic structure of the bi- 
ased BLG using the full tight-binding Hamiltonian given 
by Eq. (P). The spectrum of Eq. ^ ior V ^ reads: 



t 



V^ 



E^^iV)^±\lei + f + ^±^ti/4+{tl + V^)ei. 

(7) 

As can be seen from Eq. ^, the V = gapless sys- 
tem turns into a semiconductor with a gap controlled by 
V. Moreover, the two bands close to zero energy are 
deformed near the corners of the BZ^ so that the mini- 
mum of Ifi'k ~(^)l '^^ longer occurs at these corners. As a 
consequence, the low doping Fermi surface is completely 
different from the V — case, with its shape controlled 
byT/. 



V Vvj 
2 ^ e, 



E^-{V)^±^T^P'± ^'"-^ 



tlv 



P , 



(9) 



which explains the Fermi ring. If, instead, we have V < 
vfp <^ tj_ , we can approximate Eq. ([7]) by 



E^ 



-{V)^±^vy4 + v^pyti, 



(10) 



which corresponds exactly to the eigenvalues of the ef- 
fective two-band Hamiltonian in Eq. ([5]). Note that no 
continuum approximation can produce the trigonal dis- 
tortion shown in Fig. [T]Ja). 

The gap between conduction and valence bands, Ag, is 
twice the minimum value of Ey^~ (V) due to electron-hole 



symnietry. and is given by, 



A 



Wnyvw+y^) 



9 \ I f^ 

2t\ 9 + 



y2 
4F 



,fi+V^ 



V<Vc 



V>Vc 



(11) 

where Vc = [18t^ " ^1 + (IS^t'' + 4)^/^]^/^ ^ 6t, the 
approximation being vahd for t±^ <C t. From Eq. (jlip it 
can be seen that for both V <^ t± and V ^ t one finds 
Ag ^ V. However, there is a region for t± < V < 6t 
where the gap shows a plateau Ag ^ t±, as depicted in 
Fig. [Ijb). The plateau ends when 1/ ~ 6i (not shown). 



B. Screening of the external field 

So far we have considered V, i.e. the electrostatic en- 
ergy difference between layers felt by a single electron, as 
a band parameter that controls the gap. However, the 
parameter V can be related with the perpendicular elec- 
tric field applied to BLG, avoiding the introduction of an 
extra free parameter in the present theory. 

Let us call E = Eiz the perpendicular electric field 
felt by electrons in BLG. The corresponding electrostatic 
energy U{z) for an electron of charge — e is related to the 
electric field as eE = dU{z)/dz, and thus V is given by 



V = U{zi) - U{z2) = eEd, 



(12) 



where zi and Z2 are the positions of layer 1 and 2, respec- 
tively, and d = zi — Z2 = 3.4 A is the inter-layer distance. 
Given the experimental conditions, the value of E can 
be calculated under a few assumptions, as detailed in the 
following. 



1. External field in real systems 

If we assume the electric field E in Eq. P^ to be due 
exclusively to the external electric field applied to BLG, 
E = Eext, all we need in order to know V is the value of 



V = eEextd. 



(13) 



The experimental realization of a biased BLG has been 
achieved in epitaxial BLG through chemical doping^ii^ 
and in back gated exfoliated BLGi^i^^ In either case the 
value of E^xt can be extracted assuming a simple parallel 
plate capacitor model. 

In the case of exfoliated BLG, devices are prepared by 
micromechanical cleavage of graphite on top of an ox- 
idized silicon wafer (300 nm of Si02), as shown in the 
left panel of Fig. [UJa). A back gate voltage Vg applied 
between the sample and the Si wafer induces charge car- 
riers due to the electric field effect, resulting in carrier 
densities Ug = pVg relatively to half-filling {ug > 
for electrons and Ug < for holes). The geometry of 
the resulting capacitor determines the coefficient /S. In 



particular, the electric field inside the oxidized layer is 
Eox = e'^g/(£si02^o)i where £3102 and Eq are the permit- 
tivities of Si02 and free space, respectively. This implies 
a gate voltage Vg = eugt/ (65102^0) : from which we ob- 
tain the coefficient /3 = £si02£o/(ei). For a Si02 thickness 
t ~ 300 nm and a dielectric constant esiO-, == 3.9 we ob- 
tain /3 = 7.2 X 10^°cm~^/V, which is in agreement with 
the values found experimentallyj ^i^^i^^ In order to con- 
trol independently the gap value and the Fermi level, in 
Ref. [13 the devices have been chemically doped by depo- 
sition of NH3 on top of the upper layer, which adsorbed 
on graphene and effectively acted as a top gate providing 
a fixed electron density uqi^ Charge conservation then 



implies a total density n in BLG given by n 
or in terms of the applied gate voltage, 

n = /3Vg +no. 



no, 



(14) 



In Fig. m^b) the charge density in BLG is shown as 
a function of Vg. The symbols are the experimental 
result obtained from Hall effect measurements ^^ and 
the line is a linear fit with Eq. P^ . The fit provides 
no, which for this particular experimental realization is 
no ~ 1.8 X lO^^cm"^, and validates the parallel plate 
capacitor model applied to the back gate, since the fitted 
f3 ~ 7.2 X 10^°cm^^/V is in excellent agreement with the 
theoretical value. Extending the parallel plate capacitor 
model to include the effect of dopants, the external field 
Eext is the result of charged surfaces placed above and 
below BLG. The accumulation or depletion layer in the Si 
wafer contributes with an electric field E'b = eng/{2erSo), 
while dopants above BLG effectively provide the second 
charged surface with electric field Et = —eno/{2erSo). A 
relative dielectric constant e^ different from unity may 
be attributed to the presence of Si02 below and vac- 
uum on top, which gives e^ ~ (^Si02 + l)/2 ~ 2.5, a 
value that can be slightly different due to adsorption 
of water molecules i^^'^^ Adding the two contributions, 
Eext = E\^ + Et, and making use of the charge conser- 
vation relation, we arrive at an electrostatic energy dif- 
ference V [Eq. (fT3|) ] that depends linearly on the BLG 
density. 



y = I 2 

."0 



e^UQd 

2er£o 



(15) 



In treating the dopants as a homogeneous charged layer 
we ignore possible lattice distortion induced by adsorbed 
molecules, as well as the electric field due to the NH3 
electric dipole, which may contribute to the gap in the 
spectrum. However, it has been shown recently^ that 
for NH3 these effects counteract, giving rise to a much 
smaller gap than other dopant molecules^ as for in- 
stance NH2 and CH3. For the biased BLG realized in 
Ref. |9|, independence of Fermi level and carrier density 
was achieved with a real top gate, which makes the par- 
allel plate capacitor model a suitable approximation in 
that case. 

In the case of epitaxial BLG, devices are grown on SiC 
by thermal decomposition of the Si-faccj^ The substrate 
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Figure 2: (Color online) (a) Biased BLG devices, (b) n vs Vg 
for the left device shown in (a): experimental data is shown as 
symbols^ the line is a linear fit with Eq. (|f4p . (c) V vs n for 
the right BLG device shown in (a): symbols are experimental 
data from Ref [a; the lines are the result of Eq. (|16p . 



is fixed (SiC), and graphene behavior develops for carbon 
layers above the buffer layer,—"— as schematically shown 
in the right panel of Fig. [Ua)j^ Due to charge trans- 
fer from substrate to film, the as-prepared BLG devices 
appear electron doped with density ria- First-principles 
calculations indicate that such doping is coming from in- 
terface states that develop between the buffer layer and 
the Si-terminated substratci^i^ [Scanning tunneling mi- 
croscopy (STM) measurements corroborate the presence 
of interface stateSi ^^i^^" — ] From the point of view of our 
theoretical approach, we may interpret these interface 
states as an effective depletion layer that provides the 
external electric field necessary to make the system a bi- 
ased BLG. In Ref. [§1 the BLG density n was varied by 
doping the system with potassium (K) on top of the up- 
per layer [see Fig. [2l^a)[, which originates an additional 
charged layer contributing to the external electric field. 
Applying the same parallel plate capacitor model as be- 
fore, we get an electrostatic energy difference that can be 
written as 



V 



2er£Q 



(16) 



Following a similar reasoning to the case of exfoliated 
graphene on top of Si02, we would write Sr ~ (esiC + 
l)/2 « 5. However, this value neglects that interface 
states (the effective bottom plate capacitor) occur above 
the SiC substrate, close to the graphene system, and thus 
Ej. ~ 1 should be more appropriate. In Fig. [DJc) we com- 
pare Eq. (|16p with experimental results for V obtained 
by fitting ARPES measurements from Ref. |a. For this 
particular biased BLG realization, the as-prepared car- 
rier density was rin « 10"'^^ cm~^. From Eq. p6p . this na 
value implies a zero 1/, i.e.. zero electric field and there- 
fore zero gap, for the bilayer density n*'' « 2 x 10^'^ cm"^. 
Experimentally, a zero gap was found around n^'^^ « 
2.3 X lO^^cm"^. Given the simplicity of the theory, it 
can be said that n}^ and n'^^P are in good agreement. 
However, the agreement is only good at y ~ 0, since the 



measured V is not a linear function of n, as Eq. (|16p im- 
plies. In what follows we analyze in detail the effect of 
screening and how it modifies Eqs. (fT5|) and (|16p. 



2. Screening eorrection 

In deriving Eqs. p3|) and p^ we assumed that the 
electric field E in the BLG region was exactly the ex- 
ternal one, Eext- There is, however, an obvious addi- 
tional contribution: the external electric field polarizes 
the BLG, inducing some charge asymmetry between the 
two graphene layers, which in turn give rise to an internal 
electric field, Eint , that screens the external one. 

To estimate Eint we can again apply a parallel plate 
capacitor model. The internal electric field due to the 
charge asymmetry between planes may thus be written 
as 



Ej, 



eAn 

2£r£0 



(17) 



where —eAn is the induced charge imbalance between 
layers, which can be estimated through the weight of the 
wave functions in each layer, 



An = ni — 712 



N,Ao 



3,l = ± k 



(I^^ALkl 



' + \^'mM\'~\^'l2.\' 



<2Mn' (18) 



where the factor 2 comes from spin degeneracy, N^ is the 
number of unit cells and Aq — a^v3/2 is the unit cell 
area, jl is a band label, and the prime sum runs over all 
occupied k's in the first BZ. The amplitudes ip-'^^ j^ and 

V^Bik' with i = 1,2, are determined by diagonalization 
of Eq. ([!]) , enabling An to be written as 



NcAo 



et: 



j.i=± 



(4 + i^LM - ^i^f - (4 + ^.^^^ 



"k. + ^'-i'^k,- 



(4 + K.-)i^l - K+y + i'l + K+)tlK- 



(19) 



jl _ 



PJ^2 , 



where ek is the SLG dispersion, /C^'^ = {V/2±El^Y with 

i?^ given by Eq. ([7])- Taking the limit iVc — > oo, it is pos- 
sible to write Eq. P^ as an energy integral weighted by 
the density of states of SLG, as described in Appendix [^ 
What is important to note is that in order to calculate 
An we must specify V , which in turn depend An through 
Eq. pT)) . Thus, a self-consistent procedure must be fol- 
lowed. In particular, for the two experimental realiza- 
tions of biased BLG discussed in Sec. IIIIB 11 the self- 
consistent equation that determines V reads: in the case 
of exfoliated BLGi^ 



V 



n An(n, V) 

no no 



e^nod 

2er£Q 



(20) 
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Figure 3: (Color online) (a)-(b) Respectively, screened electric 
field and charge imbalance vs Eext at half-filling; the insets 
show the effect of changing n at fixed Eext = 0.3 V/nm, sig- 
naled by the (blue) dot in main panels, (c) V vs n for the 
BLG device shown in the right panel of Fig. [2ja): symbols are 
experimental data from Ref. |8|; lines are the result of Eq. (I21|) 
for Er = 1; the effect of changing e^ = 1 — 5 is shown in the 
inset, (d) Gap vs n for the BLG device shown in the left 
panel of Fig. [^a) with t± ~ 0.22 eV and £r = 1; the left inset 



compares the no = 5.4 x 10 cm result for 



1 (green 



dashed-dotted) with Er = 2 (blue full line); the right inset 
shows the no = 5.4 x 10^^ cm"'^ result for the screened V 
given by Eq. (|20|l (dashed-dotted line) and for the unscreened 
V given by Eq. psp (full line). We used as in-plane hopping 
t~3eV. 



in the case of epitaxial BLG;^ 



V 



n 



An{n,V) 



2erSo 



(21) 



The self-consistent electric field E = Eext + Eint at the 
BLG region, with Eint given by Eq. ([T7| for e^ = 1, is 
shown at half-filling as a function of E^xt in Fig. ^a). 
The screened E is approximately a linear function of 
Eg^t, with a constant of proportionality that depends 
on the specific value of t±. Increasing t± leads to an 
increased screening, which can be understood as due to 
an increased charge imbalance between layers, as shown 
in Fig. [3ljb) . The highly non- linear effect of inducing a 
finite carrier density (n ^ 0) can be seen in the insets of 
Fig.[3i;a) andEl^b), for t±_ = O.lt and E^xt = 0.3 V/nm. 

As a validation test to the present self-consistent treat- 
ment, we compare Eq. (j2ip with experimental results for 
V obtained by fitting ARPES measurements from Ref. |^, 
as mentioned in Sec. IIIIB II The result is shown in 
Fig. [^c) . Clearly, the self-consistent V given by Eq. ((2T|) 
for e^ = 1 is a much better approximation than the un- 
screened result of Eq. (|16p [see Fig. [D^c)]. The best fit 
is obtained for Er ^ 1 — 2, as can be seen in the inset 



of Fig. ^c). The value e^ w (esiC + l)/2 « 5 is too 
high, possibly because the bottom capacitor plate is, in- 
deed, due to interface stateS)^ and therefore is not buried 
inside the SiC substratCi ^^i^^i^^ Note, however, that the 
dielectric constant Sr may effectively be tuned externally, 
as recently shown in SLG by adding a water overlayer in 
ultra-high vacuum;^ In Fig. ^d) we show the gap Ag 
as a function of carrier density n for the biased BLG 
device shown in the left panel of Fig. HJa), with realis- 
tic values of chemical doping uqi^ The gap is given by 
Eq. mi), with ix - 0.22 eVi^ and V obtained by solv- 
ing self-consistently Eq. (PO)) for Er = I- Note that for 
Eext = we always have Eint — (the charge imbalance 
must be externally induced), and therefore we also have 
F = and Ag = 0. For this particular biased BLG device 



for n — 2rio, which 
n shown in Fig. ^d). 



the present model predicts E^^ 
explains the asymmetry for Ag 

The most important characteristic of such devices, 
from the point of view of applications, is the maximum 
size of the gap which could be induced. The maximum 
Ag occurs when Vg reaches its maximum, which occurs 
just before the breakdown of Si02. The breakdown field 
for Si02 is > IV/nm, meaning that Vg values as high as 
300 V are possible for the device shown in the left panel 
of Fig.[li;a). From Eq. HID we see that Vg ~ ±300 V im- 
plies n~ no ~ ±22 x 10^^ cm~^, and therefore Fig. [Sjd) 
nearly spans the interval of possible densities. It is ap- 
parent, specially for no = 5.4 x 10^^ cm~^, that when the 
maximum allowed densities are reached the gap seems to 
be approaching a saturation limit. This saturation is eas- 
ily identified with the plateau shown in Fig. [TJb) for Ag 
vs V, occurring for V >t±. We may then conclude that 
such devices enable the entire range of allowed gaps (up 
to tj_) to be accessed — as has been shown in very recent 
experiments.— lii The effect of using a different dielectric 
constant (e^ = 2) is shown as a full line in the left inset 
of Fig. [3t^d) , and the result for the unscreened case in the 
right inset, both for no = 5.4 x lO^^cm"^. The former 
makes the gap slightly smaller, and the latter slightly 
larger, but the main conclusions remain. 



3. Screening in continuum models 

The self consistent Hartree approach considered in the 
previous section has been applied to the full tight-binding 
Hamiltonian given in Eq. ([T]). Here we compare the re- 
sults for the potential difference V and gap Ag when 
the screening correction is used within the continuum 
approximation, either for the 4-band model of Eq. ([5]) 
or for the 2-band model of Eq. ^ . This self consistent 
Hartree approach in the continuum has been followed in 
Refs. [nil 

In the case of the 4-band model, An is still given 
by Eq. p9p with the substitutions ek —> vpp and 

^ '^ ^^ ^ ^ ■ y^ J^^ dpp, where the 



^7,i=d 



Nr^AQ ^],l = ± /_^ 



7r/i- 



prime on the right hand summation means sum over to- 



tal or partially occupied bands. Depending on the band 
in question and the value of the Fermi energy E^, the 
limits of integration 'Aiepi,p2 = {0,p , A}, where 



yr^ /(J dpp = ^— <=> A = hyJir/AQ. As regards the 



v^p^ = \JEl + V^/A ± ^El{V^ + 4) - tiFV4, 

(22) 

and A is a BZ cutoff that can be chosen such that 
471- fA _ _ " 

gap Ag , in the 4-band model it is still given by Eq. ([TT|) . 
For the 2-band model case, the charge imbalance can 
be written as an integral in momentum space of the func- 
tion \^Bi? - I0.12P = ±F/(y2 + 44pVii)i/2^ ^here 
$ — {4>bit4'A2) is the two component wave function ob- 
tained by diagonalizing Eq. ^. The ± signs stand for 
the contribution of valence and conduction bands, respec- 
tively. In particular, at half-filling the charge imbalance 
is given by 



A' 



"1/2 



t±v 



ln{2tJ\V\ + ^4tl/V^ + l), (23) 



where we have included a factor of 4 to account for both 
spin and valley degeneracies. The BZ cutoff A has been 
chosen such that up A = t^.— Since in the 2-band model 
it is assumed that V <^ t± holds we can write Aj^ /2 ~ 
-ti_V/{2TTvlh'^)\n{4tj_/\V\), which, from Eq. ^, leads 
to the logarithmic divergence of the screening ratio at 
small external electric field, E^xt/E ~ —InE, as men- 
tioned in Ref. [ij. For a general filling n the charge im- 
balance reads 



Ar 



t±v 

2iTvlh? 



In 



2i^ 




Vphf^TT^J 



4t1 



(24) 



where the charge density is given in terms of the Fermi 
wave vector as n = ±p|/(7r/i^). Inserting Eq. p4| into 
Eq. (POJ) or (|2ip we get the expression for V in the 2-band 
approximation, which is exactly the gap in the 2-band 
model, Ag = \V\. 

In Fig. mja) the obtained electrostatic energy differ- 
ence between planes V is shown for the three different 
approaches discussed above. The full (black) lines stand 
for the full tight-binding result, with V given by Eq. (PD)) 
and the charge imbalance An by Eq. P^ . The result ob- 
tained in the 4-band approximation is shown as dashed 
(red) lines. It can hardly be distinguished from the full 
tight-binding result, even when the chemical doping ng 
is as high as 5.4 x lO^^cm"^ (see figure caption). In fact, 
the only prerequisite for the continuum 4-band approxi- 
mation [Eq. (0)] to hold is that |_Ef| ^ t, which is always 
realized for the available BLG devices. As regards the 
2-band approximation model, we show as dotted (blue) 
lines the self-consistent result for V, obtained fro An as 
in Eq. ([M)) . Clearly, it is only when both the bilayer den- 
sity n and the chemical doping no are small enough for 
the relation \E-p\, V ^ t± to hold that the 2-band model 
is a good approximation (see inset). The same conclu- 
sions apply to the behavior of the gap Ag as a function 
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Figure 4: (Color online) (a) Screened V ys n for the 
BLG system shown in the left panel of Fig. [2]^a) com- 
puted within three different approaches (see text): full tight- 
binding (TB), 4-band approximation, and 2-band approxi- 
mation. Three different chemical dopings have been con- 
sidered, no = {0,1.8,5.4} X lO'^^cm"^. The inset shows a 
zoom around 1/ = for no = {0,1.8} x 10^^ cm"^ (b)- 
(c) Screened gap vs n obtained using V shown in (a), respec- 
tively for no = {0, 1.8, 5.4} x 10^^ cm~^ Parameters; i ~ 3eV, 
t± ~0.22eV, and £r = 1. 



of carrier density n, which is shown in panels |31[b)-(d) 
for no = {0, 1.8,5.4} x lO^^cm"^, respectively. The fail- 
ure of the 2-band model in the presence of interactions 
was also observed in Hartree calculations of the electron 
compressibilityi^i 



4- Electron-hole asymmetry 



As we have seen in Secs. lIIIB II and [1116 21 the two bi- 
ased BLG devices shown in Fig. [DJa) have zero gap when 
the carrier density is twice the system's chemical dop- 
ing. The closing of the gap at a finite density induces an 
electron-hole asymmetric behavior in the system, where 
obvious examples are the gap Ag and the electrostatic en- 
ergy difference between layers V, as shown in Figs, ^d) 
and|31[a). An experimental confirmation for this electron- 
hole asymmetric behavior comes from measurements of 
the cyclotron mass in the biased BLG device shown in 
the left panel of Fig. [lKa)iS (discussed in more detail in 
Sec. IIV A| . However, real electron-hole asymmetry can 
also be present in BLG due to extra hopping terms, as 
mentioned in Sec. [ill Here we study how Ag and V are 
effected by the electron- hole symmetry breaking terms t', 
74, and A, taking into account the screening correction. 

Inclusion of in-plane second-NN hopping t' leads to 
a generalized version of Eq. ^, which can be written 
as i^k.t' = H\^ — {el^f /t — 3t')l, where H\^ is given by 
Eq. (H]), ek is the SLG dispersion, and 1 is the 4x4 iden- 
tity matrix. The generalized the BLG dispersion, either 
biased or unbiased, is given by the t' ~ result added by 
—e^t'/t + 3t', which clearly breaks electron-hole symme- 
try. Note that a finite t' has no influence on the wavefunc- 
tions' amplitude. Therefore, the integrand in Eq. ([T5| - 
the definition of the charge carrier imbalance between 
layers An - is independent of t'. We have found numeri- 
cally, using a 4-band continuum model, that neither the 
screened V nor the gap Ag are affected by t', although 
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Figure 5: (Color online) (a)-(b) Band structure around K 
for the biased BLG with i' = O.li and 74 = O.lt, respec- 
tively, for V = t± = O.li. Dashed lines: t' = 74 = 0. (c)- 
(d) Respectively, V vs n and Ag vs n for the BLG device 
shown in the left panel of Fig. [2][a), modeled with a finite 
74. Parameters: t ~ 3eV, t± = O.lt, 74 = O.lt, Sr = 1, and 
no = {0, 1.8} X 10^^ cm"^. Dashed lines: t' = 74 = 0. 



the gap becomes indirect for finite t'. This means that 
the structm-e of occupied k's is insensitive to t' , and thus 
An in Eq. (fT8| is fully t' independent, at least as long as 
Ef <C t. Even though the presence of t' can lead to the 
suppression of the Mexican hat in the valence band, this 
only happens for \V\ < t']_t' ^ 10~'^t. For such a small 
\V\ value the Mexican hat plays an irrelevant role. The 
band structure around the K point for t' = O.lt (solid 
line) and t' = (dashed line) can be seen in Fig. [S{a) for 
typical parameter values. 

Now we turn to the effect of the inter-layer second- 
NN hopping 74. The generalized version of Eq. (|4]) for 
finite 74, which we call i^k,74, can be obtained by replac- 
ing the null entries {A1,A2) and {B1,B2) by 74sJ^ and 
{A2,A1) and {B2,B1) by 74Sk- The associated eigen- 
problem admits an analytic treatment at low energies 
and small biases vpp, V -^ t±r^ but as has been seen 
previously V ^ t± is possible in real systems. Therefore, 
we analyze the problem numerically using a 4-band con- 
tinuum approximation. The matrix Hamiltonian -ffk.74 
may then be written as Hk.^^ = AI^ HK.^i^^J^ near the 
K points, with M = diag[l, e*'^p, e"*'^p, 1], and HK,'yi 
obtained from Eq. ^ with iy9p = and the null en- 
tries (Al, A2), (Bl, B2), {A2, Al), and (B2, Bl) replaced 
by —V4P, where V4 = j4ah'^^^/3/2 < lO^ms"^. The 
canonical transformation defined by M clearly shows that 
the problem still has cylindrical symmetry in the con- 
tinuum approximation. Around the K' points we have 
HK'.-yi = MHK,-f4M^ . The obtained band structure for 
74 = O.lt (solid lines) and 74 = (dashed lines) is shown 
in Fig.[5jb) for typical parameter values. Note that, even 
though the gap becomes indirect for 74 7^ 0, we still have 
Ep=o = {±V/2,±y/t\ + V^ /i} as in the 74 = case. 
The screened electrostatic energy difference between lay- 



ers V for the biased BLG device shown in the left panel 
of Fig. mja) is shown as a function of the carrier den- 
sity in Fig. El^c). The result for V has been obtained 
by solving Eq. (|20p with carrier imbalance An given by 
the continuum version of Eq. (|18p. with wavefunctions 
obtained numerically through HK,~fi for 74 = O.li (see 
figure caption for other parameter values). The corre- 
sponding screened gap Ag is shown in panel [SJd) . The 
74 = result is also shown as a dashed line for both V 
and Ag . The effect of 74 may clearly be considered small, 
even for such a large value as 74 ~ 0.3 eV. However, elec- 
tronic properties which are particularly sensitive to the 
changes of the Fermi surface (like, for instance, the cy- 
clotron mass), may, in principle, be measurably affected 
by 74. We will come back to this point in Sec. lIV Al 

As regards the on-site energy A, since it is smaller than 
74 (see Sec. In]) we consider their simultaneous effect. The 
additional term in the Hamiltonian adds to the matrix 
^k,74 the contribution diag[A, 0, 0, A], and therefore the 
4-band continuum approximation for finite 74 and A may 
be written as HK,'y4,A = Hk.j^ + diag[A,0, 0, A], where 
we use the same transformation M introduced above. 
Similarly to 74, the effect of A is negligible in both V 
and An. 



C. DOS and LDOS 

Insight into the electronic properties of biased (and 
unbiased) BLG can also be achieved by studying the 
density of states (DOS) and the local DOS (LDOS) of 
the system. In particular, the LDOS can be accessed 
through scanning tunneling microscopy/spectroscopy 
measurements;^ providing a useful way to validate theo- 
retical models. On the other hand, the knowledge of the 
DOS turns out to be very useful for practical purposes, 
as it provides a way to relate the Fermi energy Ep and 

the carrier density n in the system; |n| = J„ ^ dE p2{E), 
where p2{E) stands for the BLG DOS. 

We have computed the analytical expression for the 
DOS of BLG, valid over the entire energy spectrum and 
for zero and finite bias. The expression is given in Ap- 
pendix |B] As regards the LDOS, the results have been 
obtained using the recursive Green's function method<2^ 
The DOS and LDOS of unbiased BLG has been obtained 
previously within the effective mass approximation in 
Ref. III. The effect of disorder on the DOS and LDOS 

of BLG, both biased and unbiased, has also been studied 
recently^2i29-3i^^ 

The DOS (full line) and LDOS (dashed and dash- 
dotted lines) for the biased BLG is shown in Fig.|6lja)-(b) 
for V = 0.05i. The asymmetry between the four sublat- 
tices is evident, in particular between sites Bl and A2, 
and Al and B2, which are equivalent in the unbiased 
system. Note that close to the gap edges the states cor- 
responding to positive energies have a larger amplitude at 
Bl sites, while those corresponding to negative energies 
have a larger amplitude at A2 sites. This behavior agrees 
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Figure 6: (Color online) (a)-(b) LDOS of BLG at Al/Bl and 
A2/B2 sites, respectively, for V = 0.05t and t± = O.lt. The 
total DOS is shown as a full line, (c) LDOS at A2 sites for 
n ~ {0.2,0.8,1.4} X 10^^ cm"^ and ix = O.lf. Full lines for 
numerical results and dashed lines for Eq. (|25|) . (d) LDOS at 
Ef vs n for BLG and SLG. 



with the observation that Bl and A2 are the low en- 
ergy active sites (the basis for the 2-band model), and it 
also reflects our choice of electrostatic energies in Eq. ^ : 
+V/2 in layer 1 and —V/2 in layer 2. The asymmetry 
between Bl and A2 sites can be understood with the 2- 
band continuum model, valid for vfP, V <^ t±. Defining 
the LDOS as pbi/A2{E) - i^ Ek \^bi/A2,u\^S{E - Sk), 
where $k = ('/'si.k, '/'A2,k) is the two component wave 
function obtained by diagonalizing Eq. (JG)), we can read- 
ily arrive at the following expressions. 



PBl/A2iE) 



1 t^ 

2V37r t^ ' 



ME) 



E±V/2 



(25) 



The asymmetric behavior is apparent, with pbi{E) di- 
verging for E -^ V/2+ and pA2iE) for E -^ -V/2~ . 
The result for pA2{E) is shown in Fig. ^c) for V ~ 
{0.87,4.23,7.87} x IQ-^f and U = O.li. Within the 
screening corrected parallel plate capacitor model dis- 
cussed in Sec. lIlIB] [Eq. (PU)) ]. these V values correspond 
to carrier densities n ~ {0.2,0.8,1.4} x lO^^cm"^, re- 
spectively, where we have used t ~ 3.1 eV, tt-q == 0, and 
£r = \- The full lines are the recursive Green's func- 
tion methodic results and dashed lines are the results of 
Eq. (P5)l . As expected, the closer to the gap edges the 
better the agreement between the two approaches. 

A strong suppression of electrical noise in BLG has 
been reported recently by Lin and Avourisj^ In devices 
made from exfoliated BLG on top of Si02, the current 
fluctuations are thought to originate from the fluctuating 
trapped charges in the oxide. Therefore, the more effec- 
tive the impurity charge screening in the system the lower 
the electrical noise. The lower noise in BLG than in SLG 
may then be attributed to the low energy finite DOS in 
the former. However, it has also been reported in Ref. l76l 
that while increasing the carrier density in SLG leads to 
lower noise, as expected due to more effective impurity 



screening, it results in higher noise in BLG. Insight into 
this behavior is achieved by analyzing the LDOS at the 
Fermi level Ep in a biased BLG, as charging the system 
through the back gate Vg leads to a finite perpendicu- 
lar electric field. In Fig. [SJe) we show the biased BLG 
LDOS at E-p for Bl and A2 sites as a function of carrier 
density n in the system. For a given n, the electrostatic 
energy difference V is evaluated self-consistently through 
Eq. (P0| . with riQ = and e^ = 1, and E-p is obtained by 
integrating over the DOS. Additionally, we use i ~ 3eV 
and t±_ = O.li. We have chosen densities in the range 
n € [0 — 2] X lO^^cm"^, which corresponds to back gate 
voltages Fg e [0 - 27] eV through Eq. dM]) , similar to the 
experimental range in Ref. l76l . The main observation to 
be made as regards the results of Fig. ^e) is that for the 
low energy active sublattices Bl and A2 the LDOS at Ep 
remains approximately constant with increasing electron 
density, as opposed to the ~ ^/n dependence found in 
SLG. This is an indication that impurity screening may 
not be increasing with carrier density in the biased BLG 
system, which may be contributing to enhance electrical 



IV. MAGNETIC FIELD EFFECTS 

In the biased BLG system, as a consequence of the 
gapped band structure discussed in Sec. IIIIl a perpen- 
dicular magnetic field is expected to induce distinct fea- 
tures in electronic properties. In this section we focus 
on the cyclotron mass (semi-classical approach) and on 
the cyclotron resonance (quantum regime) comparing the 
theory with experimental results. 



A. Cyclotron mass 

In the semi-classical approximation the cyclotron ef- 
fective mass 771c is given by 



h^ dA{E) 



2ti dE 



E — £/p 



(26) 



where A{E) is the fc-space area enclosed by the orbit of 
energy E, and the derivative is evaluated at the Fermi en- 
ergy Ep jZii^ It can be accessed experimentally through 
the Shubnikov-de Haas effect, providing a direct probe 
to the Fermi surface. In the case of exfoliated graphene, 
either SLG or (un)biased BLG, the Fermi energy can be 
varied by tuning the back gate voltage, and therefore a 
significant portion of the whole band structure may be 
unveiled. In particular for the biased BLG, the presence 
of a finite gap can be checked and the model developed 
in Sec, mil tested. 
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1. Comparison with experiment 

General expressions for TOc obtained for the full tight- 
binding bands in Eq. ([7]), valid for the relevant param- 
eter range V < t^ <^ t and restricted to E-p < t, 
are given in Appendix O In Fig. [Tl^a) we compare the 
theory results for the cyclotron mass with experimental 
measurementsiS on the biased BLG system shown in the 
left panel of Fig. Efa). We have only considered nic as- 
sociated with low energy bands E^~ [see Eq. ([T])], since 



0.09 



E^~^ are inactive for the experimentally available carrier 



densities. The dashed lines stand for the unscreened re- 
sult, where V is given by Eq. P^ . and the solid lines 
are the screened result, with V given by Eq. ^U^. The 
inter-layer coupling t±_ has been taken as an adjustable 
parameter, keeping all other fixed: t ~ 3eV, Sr = I, and 
no = 1.8 X 10^^ cm^^. The value of t± could then be cho- 
sen so that theory and experiment gave the same rUc for 
n = 2no « 3.6xl0^^cm~^. As discussed in Sec. fill B 21 at 
this particular density the gap closes, meaning that the 
theoretical value becomes independent of the screening 
assumptions. We found t± « 0.22 eV, in good agreement 
with values found in the literature. The theoretical de- 
pendence mc{n) agrees well with the experimental data 
for the case of electron doping. Also, as seen in Fig.[71[a), 
the screened result provides a somewhat better fit than 
the unscreened model, especially at low electron densi- 
ties. This fact, along with the good agreement found 
for the electrostatic energy difference data of Ref. |^ [see 
Fig. [3Uc)[; allows us to conclude that for doping of the 
same sign from both sides of bilayer graphene, the gap 
is well described by the screened approach. In the hole 
doping region in Fig. [7]Ja), the Hartree approach under- 
estimates the value of rric whereas the simple unscreened 
result overestimates it. This can be attributed to the 
fact that the Hartree theory used here is reliable only if 
the gap is small compared to t±. In the experimental 
realization of Ref. [l3; no > and, therefore, the theory 
works well for a wide range of electron doping n > 0, 
whereas even a modest overall hole doping n < cor- 
responds to a significant electrostatic difference between 
the two graphene layers. In this case, the unscreened 
theory overestimates the gap whereas the Hartree calcu- 
lation underestimates it. 

In Fig. Eljb) we compare our best fit to the cyclotron 
mass (full line) with results obtained for different param- 
eter values. The dashed-dotted lines stand for m^ ob- 
tained with Er = 2 in Eq. (P0|) . As can be seen clearly, the 
n > result is not substantially affected, while for n < 
the theory description of nic worsens. This is due to the 
reduction of the gap when e^ is increased [see left inset 
in Fig.[3Id)[. The dashed lines in Fig.[71[b) are obtained 
with no — 0, where the zero gap occurs at the neutrality 
point. The dotted lines are the result for E^xt = = V, 
i.e., zero gap at every density value. Note that these two 
results, no = and V = 0, show an electron-hole sym- 
metric nic, contradicting the experimental result. It may 
then be said that the electron-hole asymmetry observed 




screened 
unscreened 

I L 

-4 4 

nfio'^cml 



n(10'^cm"^) 



Figure 7: (Color online) Cyclotron mass vs n, normalized to 
the free electron mass, rria. (a) Solid lines are the result of 
the self-consistent procedure and the dashed lines correspond 
to the unscreened case; t ~ 3eV, t± ~ 0.22 eV, Sr ~ 1, and 
no = 1.8 X lO^^cm"^. Circles are experimental data from 
Ref. [l^. (b) The screened result in (a) is coinpared with the 
result for Sr = 2, the case without chemical doping (no — 0), 
and the case where the external field is zero {V = 0). 



in rric is a clear indication of the presence of a finite gap 
in the spectrum. It will be shown in Sec. IIV A3l that. if 
we ignore the gap, this electron-hole asymmetry cannot 
be described by taken into account i', 74 or A. 



2. Cyclotron mass in continuum models 

Here we compare our results for the cyclotron mass, 
which has been obtained with expressions shown in Ap- 
pendix [C] with the results of continuum models. 

Within the 4-band continuum model given by Eq. ([5|) , 
where the dispersion is just the full tight-binding result 
[Eq. ([7])[ with the substitution ek — > vfP, we can easily 
derive the following analytical expression for mc. 



Ev 



nic 



V^+tl 



2^Ej{V^+ti)-tlVyA 



(27) 



In Fig.[5Ka) the dashed line is the result of Eq. (P7|) . where 
V has been computed self-consistently using Eq. ^U\i 
and the 4-band continuum approximation discussed in 
Sec. nil B 31 As expected, the agreement with the 
full tight-binding result (shown as a full line) is ex- 
cellent for the considered densities. Note that there 
is an extra solution given by nicWp = Ef[1 — {V'^ + 
t\)I^AEl{V^ + t\) - t\V^], valid when |^f| < V/2 or 
|£'f| > •\/y^/4 + t\, which corresponds to the extra or- 
bit appearing when Ey falls in the Mexican-hat region, or 
above the bottom of high energy bands. We can estimate 
the densities for which these two regions start playing a 
role: using V ^ O.lt^ ^ O.Olf in the Mexican hat region 
(vahd for no < 2 x 10^^ cm^^) we get n < lO^^cm^^; 
setting V ^ t^ ^ {).\t around the bottom of high energy 
bands we get n > lO-'^'^cm"^. These two density values 
are outside the range of experimentally realized densities 
[see Fig.[7Ka)[. 
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Figure 8: (Color online) Cyclotron mass vs n, normalized to 
the free electron mass, rric. (a) Comparison between full tight- 
binding (TB) and 4-band approximation for t± ~ 0.22 eV and 
no = 1.8 X 10^^ cm"^ (b)-(c) Effect of finite t' and 74 for no = 
and no = 1.8x lO^'^ cm~'^, respectively: dotted line is for t' ~ 
0.3 eV and t± ~ 0.22 eV; dashed line is for 74 ~ 0.12 eV and 
t± ~ 0.19 eV; full thin line is for i' = 74 = and ix = 0.22 eV. 
(d) Effect of A for 74 ~ 0.12 eV and tj_ ~ 0.19 eV: full line 
for A ~ 0.03 eV; dotted-dashed line for A ~ -0.03 eV; dashed 
line for A = 0. Circles are experimental data from Ref. [lOl . 
We have used t ~ 3 eV and Sr = 1. 



distortion of the energy bands shown in FiglSJb). This 
attenuation can be understood as the result of fixing the 
carrier density n and not the Fermi energy Ep : changing 
74 (or t') for a given n leads to a different Ep, and the 
new Ep is such that it counteracts the expected effect of 
74 (or t') in 771c. Fig. [5]Jc) shows the same as [Sl^b) for 
TT-o = 1.8 X lO^'^ cm~^. The effect of the on-site energy A 
is shown in Fig.IiJd) for fixed 74 ~ 0.12 eV, ij_ ~ 0.19 eV 
and no = 1.8 X 10^^ cm~^. The result for A = (dashed 
line) is shown along with the result for A ~ 0.03 eV (full 
line) and A ~ — O.OSeV (dotted-dashed line). It is clear 
that the effect of t', 74 and A on the cyclotron mass can 
be neglected. 



B. Cyclotron resonance 

The effect of a perpendicular magnetic field can be 
studied within the continuum approximation through 
minimal coupling p — > p — eAX The case of biased BLG 
has been studied both within the 4-band [Eq. dSJl and 
2-band [Eq. ^] continuum models in Refs. [niil2,.^.^. 
Here we use the same approach to study the cyclotron 
resonance (i.e. the Landau level transition energies) with 
the extra ingredient that the parameter V depends on the 
filling factor, as discussed in Sec. IIIIBI 

In the 4-band model standard manipulationsiiiii^i^ 
lead to the unbiased BLG Landau level spectrum 



3. Effect of electron-hole asymmetry 



E: 



±± 



In Sec. nil B 41 the effect of electron-hole symmetry 
breaking parameters - namely, i', 74, and A - has been 
studied regarding the self-consistent description of the 
gap. Here we extend the analysis to the cyclotron mass, 
restricting ourselves to the biased BLG device shown in 
the left panel of Fig. ^a.). Results have been obtained 
within the 4-band model. As all cases have cylindrical 
symmetry around K and K' , the cyclotron mass may be 
written as nic = pp / {dEp / dpp). 

In Fig. [5]Jb) we show the nic result for finite t' (dot- 
ted red line) and finite 74 (dashed blue line), keeping 
TT-o = (absence of electron- hole asymmetry due to chem- 
ical doping). The thin full line is the result obtained for 
i' = 74 = in Sec. IIV All and circles are experimental 
data from Ref. llO- The n > region, where the smaller 
gaps are realized experimentally, is still well fitted if we 
choose i_L ~ 0.22 eV with t' ~ 0.3 eV or tx_ ~ 0.19 eV with 
74 = 0.12 eV (we use t ~ 3eV). However, it is clear that 
neither of these results can account for the electron-hole 
asymmetry observed experimentally. In fact, a closer 
look reveals that the nic for finite t' have the opposite 
trend, being smaller than the t' = result for n < 
and larger for n > 0, as would be expected by inspec- 
tion of the energy bands in Fig. [SJa) . Such an opposite 
trend should also be seen for finite 74, although the ef- 
fect is not as large as expected from the considerable 



r 



±W(l + 2n)4- + 



f2 

2 



± 



^i^^+tir/A + n^Hl, 

(28) 

where 7 — \/2vph/lB, with Ib — \/h/\e\B for the 
magnetic length. Non-zero {n > 1) eigenenergies are 
fourfold degenerate due to valley and spin degeneracy, 
while zero energy Landau levels have eightfold degen- 
eracy, since there are two zero energy Landau states 
in = —1,0) per valley per spin. The 2-band model re- 
sult E^ sa ±7^i]|"'"-\/n(n + 1) is easily recovered from 
Eq. ([28| for 7 ^ <j^, being valid for magnetic fields up 
to B « 1 T^ 

The Landau level transition energies in BLG have 
been recently obtained through cyclotron resonance 
measurements.— The data was found to deviate from 
what would be expected through Eq. (pS)) . especially for 
larger filling factors. It should be noted, however, that in 
order to keep a constant fil ling factor and vary the mag- 
netic field, as done in Ref. |82, the back gate voltage Vg 
has to be tuned to compensate for the variation of Lan- 
dau level degeneracy. As we have seen previously, tuning 
Vg is equivalent to change V - the electrostatic energy 



difference between layers - which means that Eq. ([281) is 
no longer valid, as recently shown within the 4-band con- 
tinuum modeli^ To have an estimate for the effect of the 
back gate voltage on the Landau level spacing we have 
computed Landau level energy differences taking into ac- 
count the variation of V with carrier density n. We have 
used the unscreened result given by Eq. (jlSp , with tiq = 
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Figure 9: (Color online) Landau level transition energies vs 
magnetic field for the given filling factors. Tlie dashed line is 
the unbiased BLG result [Eq. (|28|) ] and the line with crosses 
is the biased BLG result (see text). We used t = 3.5 eV and 
tj_ — Q.lf. Filled symbols are experimental data from Ref . |83: 
circles for electrons and squares for holes. 



and £,. = 1. Within this approximation we can easily 
write V in terms of the filhng factor v and magnetic field 
BasV = iyBe'^d/{2eo(j)o) « 7.4 x lO'^iyB, with B in 
Tesla in the last step. Thus, for fixed filling factor, V 
varies linearly with B. Note that the comparison between 
this unscreened treatment of the biased BLG and the 
unbiased result in Eq. (|28|) should give lower and upper 
limits for the effect of the perpendicular external field in 
the cyclotron frequency. In Fig. ^ we show the obtained 
Landau transition energies vs magnetic field for the given 
filling factors. The dashed lines represent the unbiased 
BLG result, as given by Eq. (|^5|) . The lines with crosses 
are the results for the unscreened biased BLG, and filled 
symbols are experimental data from Ref. |82| : circles for 
ly > and squares for ly < Q. We have used t = 3.5 eV and 
t± = O.lt, consistent with Ref. |82. As can be seen from 
Fig. [HI the back gate induced electric field gives rise to siz- 
able effects already for magnetic fields and filling factors 
realized in experiments. Except at j/ = ±8, the result of 
Eq. (|^5|) for the unbiased BLG and the unscreened biased 
BLG result effectively provide upper and lower limits to 
the experimental data. The observed electron-hole asym- 
metry could then be interpreted as due to an asymmetry 
in V vs n: larger V, and therefore larger gap, for n < 0; 
smaller V and gap for n > 0, which would make the re- 
sult more close to the unbiased case. It should be noted 
that in such a case we would expect the neutrality point 
to occur for Vg < 0, as is the case of the NII3 doped BLG 
studied before. For the system reported in Ref. [82, how- 
ever, the opposite seems to be happening, as indicated 
by the Hall resistivity. A neutrality point for Vg > is, 
in fact, the more usual effect of H2O molecules adsorbed 
on graphene samples.— As a final remark regarding the 
results presented in Fig. [HI we note that the experimental 
data trend, which makes Eq. (pS)) a poor fit at \iy\ > 8, 
is still not accounted for in the biased BLG result. An 



alternative approach is the inclusion of the screening cor- 
rection, which should go beyond Eq. ([T^ including the 
magnetic field effect. It has been reported recently that 
Dirac liquid renormalization may also be contributing to 
the observed trend 1^ 



V. CONCLUSIONS 

We have studied the electronic behavior of bilayer 
graphene in the presence of a perpendicular electric field 
- biased bilayer ~ using the minimal tight-binding model 
that describes the system. The effect of the perpendic- 
ular electric field has been included through a parallel 
plate capacitor model, with screening correction at the 
Hartree level. We have compared the full tight-binding 
description with its 4-band and 2-band continuum ap- 
proximations, and found that the 4-band model is always 
a suitable approximation for the conditions realized in ex- 
periments. Also, we have studied the effect of electron- 
hole asymmetry terms and found that they have only a 
small effect on the electronic properties addressed here. 
The model has been applied to real biased bilayer de- 
vices, either made out of SiG^ or exfoliated graphenei^iiS 
The good agreement with experimental results - namely, 
for the electrostatic energy difference between layers ob- 
tained through ARPES^ and for the Shubnikov-de Haas 
cyclotron massif - clearly indicates that the model is 
capturing the key ingredients, and that a finite gap is 
effectively being controlled externally. Analysis of recent 
experimental results regarding the electrical noise^^ and 
cyclotron resonance^ further suggests that the model 
can be seen as a good starting point to understaird the 
electronic properties of graphene bilayer. 
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Appendix A: Asymmetry between layers 



In order to write Eq. p9p as an energy integral, we 
start by introducing the SLG density of states per spin 
per unit cell defined for the conduction band as 



p(') = ^T.si^~t\ 



Sk 



(Al) 



with Sk as in Eq. (|4]). The momentum sum in Eq. (jAip 
can be written as an integral by letting N^. -^ 00. The in- 
tegral can be performed and written in terms of complete 
elliptic integrals of the first kindi^ 
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With the definition of p(e) in Eq. (jAip the charge im- 
balance between layers in Eq. P^ can be written as 
An = Ani/2 + An. where the charge imbalance at half- 
filling Anj^ /2 is given by 



Am 



/2 



—y 



3t 



dep{e)I-i{e), 



(A2) 



and the fluctuation An with respect to the half-filled case 
is given by 



An 



^O \-E,=±/e?dep(e)Z_,(e), n<0 



, (A3) 



where n is the carrier density with respect to half-filling. 
The integral kernel in Eqs. (jA2p and (jA3p reads 

(A4) 



where /Ci'(e) = [V/2 ± Ei\t)f , with £:J'(e) given by 
Eq. ([7]) with the substitution ci^^ e. The limits of inte- 
gration in Eq. (|A3|) depend on the band label I and Ey as 
follows: with I = — we have ei = e~ and €2 = e"*" for £^| < 
y^/4, while for £'| > y2/4 we have ei = and £2 = 6"'"; 
with ^ = + we only have contribution for £"1 > i^ + V^^/4, 
and the limits are ei = and €2 = e~. We use the nota- 



tion e 



± 



[El 



1/2/4 ± ^El{V^+i\)-i\V^lA]h. 



Appendix B: Bilayer DOS 

The DOS per unit cell of BLG, either biased or unbi- 
ased, is defined as 



P2[ 



(^) = W E['5(^ - ^k ") + 5{E i?±+)], (Bl) 



where E^ is given by Eq. ([7]). Equation (jBip can 
be written as a sum of two contributions, P2{E) = 
X]/=± PiiE), where the label / = ± stands for contribu- 
tions coming from bands E^^ . The analytical expressions 
for each contribution are 



P2iE) 



and 



ib—iE) 



X-{E) 



-K 



y/Flx-(E)/t]\F[x-{E)/t]J^ 



( ix-{E)/t \ 
\Flx-{E)/t]) 



^ ^^WFlx+(E)/t][Fix+{E)/t\ 



t^n' 



^—{E) 



X-jE) 



-K 



( F[x-{E)/t] \ 
[ 4x-{E)/t J 



Vix-{E)/t \ ix-(E)/t 
+ 






(E)/t] 



(E)/t 



Ag/2 < \E\ <V/2Aa<t^ 

V 
E+-{t) < \E\ <V/2Aa>t'^ 

Ag/2< \E\ < E+- (t) A a < t"^ 
Ag/2< \E\ <E+-{t)Aa>e 

' E+-{t) < \E\ < E+-{M) Aa<e 

V 
Ag/2 < 1^1 < S+-(3i) A t^ < a < 9t2 



pt{E) = 



i'+-iE)- 



X-jE) 



-K 



(^ 



4x'(g)/t 



^tl+Vyi<\E\<E++{t) 



t^ 



_ y/F[x-(E)/t] \F[x-(E)it]) ' 



(B2) 



(B3) 



with ?/'(£') given by 



i^^\E) = 



V4/4 + [ij + V^)x'{EY^x^EY + tl/2 + V^/A ± ^4/4 + [ij + V^)x'{E) 



I 



(B4) 



and x^{E) as in the right-hand side of Eq. (P^ with K(m) for the complete elliptic integral of the first kind, 
Ef -^ E. We use F{x) = (1 + x)'^ - (a:^ - 1)2/4 and 
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and E^^{x) is given by Eq. ^ with the substitution 
Ek ^ a; and a = {V^/4 + tlV^/2)/{V^ + tl). 



Appendix C: Cyclotron mass 

Based on the full tight-binding band structure E^ 
given in Eq. ([7|), it is possible to derive general expres- 
sions for the cyclotron mass in Eq. ([^5]) . The key ob- 
servation is that the area of a closed orbit at the Fermi 

El 
Ak, where 
k 

the prime means summation over all k's inside the or- 
bit, and Ak = {2tt)'^ / (NcAq) is the area per fc-point in 



the first BZ. The cyclotron mass may then be written as 

mc oc dE, E^ It '^^ Sk SiE " ^r)0(e^ - £k), where 
ek is the SLG dispersion and e^ is given by Eq. ([7]). The 
integration limits £i and the choice between the two pos- 
sibilities e^ depend on the particular band and on the 
position of the Fermi level. Skipping the details of the 
derivation, what is worth noting is that, due to the sum 
of delta functions in the previous expression for ttIc, the 
result has a mathematical structure similar to the de- 
rived expressions for the DOS of BLG (see Appendix iBj) . 
The cyclotron effective mass of the biased BLG for the 
relevant parameter range V < t± <^ t and \Ep\ < i is 
then given by 



mciEp) 



Aot'^TT 



-4' 



(Ef) 



X-(Ef) 



+ iEF] 

-{Ef] 



V-F[x-(BF)/i] 



X+JEf) 



-K 



y/F[x+iEF)/t] 
^F[x-{EF)/t] \ 



^( i-?x-(t)Al )' \/2<\Ef\<V/2 

( 4x+(gF)/t A 
\F[x-^[EF)/t\) 



4x-(gF)/t 
F[x-{EF)/t] 

4x+(gF)/t 
F[x+(EF)/t\ 

4x~iEF)/t 
F[x-{EF)/t] 



r 



Ag/2 <\EF\<t 



(CI) 



V^+VyI<\EF\<t 
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